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ABSTRACT 

$ZL(' Context. The Hall drift, namely, the transport of magnetic flux by the moving electrons giving rise to the electrical 

current, may be the dominant effect causing the evolution of the magnetic field in the solid crust of neutron stars. It is 
a nonlinear process that, despite a number of theoretical efforts, is still not fully understood. 

Aims. Through mostly analytic arguments and solutions, we intend to help understand this highly nonlinear process 



o 



^ ' Methods. We use the Hall induction equation in axial symmetry to obtain some general properties of nonevolving fields, 

' as well as analyzing the evolution of purely toroidal fields, their poloidal perturbations, and current-free, purely poloidal 

fields. We also analyze energy conservation in Hall instabilities and write down a variational principle for Hall equilibria. 
, Results. We show that the evolution of any toroidal magnetic field can be described by Burgers' equation, as previously 

^ ' found by Vainshtein and collaborators in a plane-parallel geometry. This evolution leads to sharp current sheets, which 

dissipate on the Hall time scale, yielding a stationary field configuration that depends on a single, suitably defined 
[ coordinate. This field, however, is unstable to poloidal perturbations, which grow as their field lines are stretched by 

, the background electron fiow, as in the instabilities found numerically by Rheinhardt and Geppert. On the other hand, 

current-free poloidal configurations are stable and could represent a long-lived crustal field supported by currents in the 
fluid stellar core. There may be additional, stable configurations, corresponding to restricted local minima or maxima 
\ of the magnetic energy. 

>— ' , Conclusions. Hall equilibria can be described by a simple variational principle. Long-lived, toroidal fields are not expected 

in neutron star crusts or other regions where Hall drift is the dominant evolution mechanism. However, other stable 
^—■^ ' configurations do exist, such as current-free poloidal fields and possibly others. 
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. 5t 1- Introduction the electrons are the only mobile charges (iJonesl 19881 : 



iGoldreich fc ReiseneggeiHi992l : [Cumming et al.l . l2004i r 
The Hall drift, namely the advection of magnetic flux by the The evolution of the magnetic field under the Hall 
current associated with it, is important in systems where effect is a nonlinear process that has so far eluded a 
the magnetic field is strong enough to make the cyclotron full theoretical understanding. Since the Hall effe c t con- 
frequency comparable or greater than the collision rate, serves the magnetic eii e rgv (lUrpin &: ShalvbkovL Il99ll : 
and bulk flow velocities are not much larger than the rel- IGoldreich fc Reiseneggeii Il992f) . it has been argued that 
ative velocity of different charge carriers, associated to the it will act mainly through the generation of steep mag- 
electric current. For some time, it has been realized that netic field gradients, on which resistive dissipation can 
this effect might play an important role in the evolution act much more quickly than it would on a smooth field, 
of neutron star magnetic fields, which under some circum- IGoldreich fc Reiseneggeii (|1992t ) gave general arguments 
stances are large enough for the above conditions to be sat- showing that this could happen through a turbulent cas- 
isfied, particularly in the solid crust of these stars, in which cade transferring energy from larger to smallcr scalcs, which 

was later supported by simulations (|BiskamD et al.l . Il999t ). 



On the other hand, I Vainshtein et al.l (|2000| ) considered an 
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pation is included, but indeed lead to rapid dissipation of 
the magnetic field. 

Other attempts at a bett er understanding have 
been numerical. Sortie authors (i Naito fc Koiimal . 119941 : 
iHollerbach fc Rtidigeii [2002112004 have computed the evo- 
lution of the first few eigenmodes of the magnetic diffusion 
equation in spherical symr netry, nonlinearly co u pled t o each 
other by the Hall effect. lUrpin fc Shalvbkovl (|l99lt ) have 
modelled the evolution of the magnetic field by solving the 
evolution equation on a grid for a purely toroidal field on 
a uniform density sphere, again finding that the magnetic 
field develops sharper gradients, which later dissipate by re- 
sistive effects. Simulating the evolution o f a poloidal field in 
a similar way (IShalvbkov fc UrDinl . [T996h . energy was found 
to be transferred to a toroidal co mponent and back , com - 
plicating the results substantially. iPons fc GeppertI (|2007l ) 
have recently simulated the evolution of both poloidal and 
toroidal components in a neutron star crust with a code 
that combines finite differences in the radial direction with 
a spherical harmonic decomposition in angle, using a re- 
alistic description of the radial profile of electron density 
and electric resistivity and the evolution of the latter. They 
confirm the transfer of energy back and forth between the 
poloidal and the toroidal component (especially from the 
former to the latter), and find regimes along the neutron 
star evolution when Hall drift or Ohmic diffusion dominate. 

In a different approach, [R heinhard t~fc GeppertI (12002 ) 
(see alsolGeppert fc Rlieinhardtii2002tiGeppert et al.]|2003 : 



iRheinhardt et all 120041 )" did a numerical stability analy- 
sis on plane-parallel field configurations that do not by 
themselves evolve under the Hall effect, finding both sta- 
ble modes (a particular case of which are the well-known 
"whistler waves") and unstable ones. Based on several ex- 
amples, they conjecture that, at least in simple geometries, 
the determining factor for the existence of unstable modes 
is a sufficiently large second spatial derivative of the unper- 
turbed magnetic field strength. 

The present paper reports an analytical study of the 
evolution of a magnetic field with a possibly realistic, ax- 
ially symmetric configuration in a solid star with non- 
uniform electron density. We first give a short re-derivation 
of the relevant form of the induction equation (§ [5]) and 
specialize it to the two independent scalar functions de- 
termining the magnetic field in axial symmetry, discussing 
some of its implications (§[21)- Then, we consider the very 
special case of a purely toroidal magnetic field (§ Hj). It 
is shown that, if only the Hall drift is active, the mag- 
netic field evolves independently on each of a family of 
nested toroidal surfaces, again being determined (for a suit- 
able cha nge of variables) by Bur gers' equation 14. ip . As 
found bv IVainshtein et all (j2000t ). this leads to the forma- 
tion of current sheets that can rapidly dissipate, so the 
field evolves on a Hall time scale to a simple configuration 
on each of these surfaces, minimizing the magnetic energy 
while conserving the flux within each toroid (§ I4.2p . Next 
(§0; we consider various aspects of small perturbations and 
Hall instabilities: We show that a stationary toroidal mag- 
netic field is generally unstable to poloidal perturbations 
(ij lS.ll — up to here, our r esults were already summarized 
in iReisenegger et "aLll2005D and interpret this through the 
plane-parallel model considered by previous authors fij |5.2p . 
we discuss the issue of energy conservation in Hall instabili- 
ties, show that current-free configurations are stable fi) l5.3p . 
and that a current-free, poloidal field supports whistler-like 



waves fii l5.4p . Then, we show that Hall equilibria can be de- 
rived from a variational principle, by requiring the magnetic 
energy to be stationary with respect to small displacements 
that do not alter the electron density, and use this to ar- 
gue that some configurations, corresponding to maxima or 
minima of the magnetic energy subject to these constraints, 
must be stable with respect to small perturbations evolving 
solely through the Hall drift (§111). Finally, we list our main 
conclusions (§|2|). 



2. Induction equation 

In order to clarify our assumptions and notation, we start 
by re-deriving the equation of motion for the magnetic 
field in an arbitrary geometry. More general derivations 
have been given, among other s . inlGoldreich fc Reiseneggeil 
(|1992D and IReisenegger et all ()2005D . 

We consider a material in which the electrons are the 
only moving particles, embedded in a perfectly rigid, neu- 
tralizing background (i. e., an idealized crystal lattice) 
against which they can occasionally scatter. The steady- 
state, local average velocity of the electrons is determined 
by the balance of the Lorentz force against the time- 
averaged momentum loss through collisions, which yields 
a generalized Ohm's law. 



-new 



a E 



V X B 



where j is the electric current density, n, — e, and v are the 
number density, charge, and average velocity of the elec- 
trons, (J is the (scalar) electrical conductivity, E and B are 
the electric and magnetic field, and c is the speed of light. 
The evolution of the magnetic field with time t is described 
by the induction equation, dS/dt — — cV x E, which, re- 
placing E from eq. (1), takes the form 



9B 

'dt 



V X V X B 



cj 



(2) 



In the infinite-conductivity limit, this equation asserts that 
the magnetic field lines can be thought of as drifting along 
with the electron flow ("Hall drift"). The term involving 
the conductivity corresponds to deviations from this ideal- 
ization due to resistive diffusion of the magnetic fleld. 

In the slow- motion limit (i. e., ignoring electromagnetic 
radiation), the current, and thus the electron velocity, are 
also related to the magnetic field by Ampere's law, j — 
—nev — cV X B/ (47r) . This yields a magnetic fleld evolution 
law with B as the only dynamical variable. 



'dt 



-V X 



Airne 



(V X B) X B + 77V X B 



(3) 



where we have introduced the magnetic diffusivity, rj = 
c^/(47rcr). The Hall term is quadratic in B, which has so 
far impeded a full theoretical understanding. 



3. Axial symmetry 

We now specialize to the case of an axisymmetric star, de- 
scribed by the standard cylindrical coordinates R, (j), and z, 
with electron density n{R, z). For simplicity, we set 77 = 
in this derivation, although below we consider the effects 



A. Reisenegger et al.: Hall drift 



3 



of resistive diffusion. The most general, axisymmetric mag- 
netic field can be decomposed into a toroidal component 



Bt = B{R,z)Vcj) 

and a poloidal component 

Bp = VA{R,z) X V0. 



(4) 



(5) 



This decomposition makes it explicit that the field depends 
only on two scalar functions, B and A^ and explicitly satis- 
fies the condition of zero divergence independently for both 
components. We choose to write it in terms of the gradient 
\J(j) = (pl R instead of the unit vector (f> in order to make 
easy use of the identity V x V(/) — 0. For future reference, 
we also write the toroidal and poloidal components of the 
electron velocity. 



Vp = 



47ren(i?,2:) 



47ren(i?,2:) 



VB X V0, 



(6) 
(7) 



where A* = - (2/R)d /dR is the Grad-Shafranov oper- 
ator (e. g.. lKulsrudl[2005l § 4.9). 

This shows that the magnetic field lines lie on the sur- 
faces A = constant, while the current or electron flow lines 
lie on surfaces B = constant. If both A and B are taken to 
be zero on the symmetry axis, then 2TrA is the poloidal flux 
enclosed by a given surface A = constant, whereas cB/2 is 
the total current e nclosed by the corresponding surface (see 
also iKulsrudI [20051 . § 4.9). 

Replacing the previous expressions into the induction 
equation, eq. ([2]), and separating toroidal and poloidal com- 
ponents, we obtain the evolution of the two scalar functions, 



= i?2 [v(xAM) xVA + Vxx BVB] ■ V(/), (8) 



dB 

'di 
dA 

'dt 



= R^x^A X VB ■ V0, 



(9) 



where we have introduced the scalar function x(^j ^) = 
c/[4:TTen{R, z)R'^]. These equations show that the toroidal 
and poloidal components of the field couple strongly to each 
other in their evolution, so the latter might become quite 
intricate in the general case. However, there are a few simple 
but interesting things to be noted. 

First, A does not evolve if the surfaces A = constant 
and B = constant coincide, i. e., one of these scalar func- 
tions can be written as a function of the other. In this 
case, the poloidal components of the magnetic field and the 
current density are parallel to each other at every point, 
jp — {c/ATT){dB/dA)'Bp. They lie in these surfaces, and 
the coefficient is also constant on each surface. Since there 
cannot be substantial currents flowing into and out of the 
star, the poloidal currents on any flcld line extending out 
of the star must vanish identically. 

This condition, jp x Bp 0, is equivalent to imposing 
that the azimuthal component of the Lorentz force is zero. 
Since neutron stars are born in a fluid state, it might be 
natural to require that they start in an MHD equilibrium, 
in which magnetic forces are balanced by fluid forces. In 
an axially symmetric configuration, the fiuid forces cannot 
have an azimuthal component, therefore naturally leading 
to this condition. 

Second, the previous condition does not ensure that B 
will not evolve, so it may not continue to be satisfied at later 



tim es as the magnetic field Hall-d rifts through the solid 
(see Reisene gger fc Thompsonll200"7l for a related discussion 
for Hall drift in a fluid star). In order to insure that B = 
constant as well, one also needs to impose that x(^*-^ + 
BdB/dA) is constant on the same surfaces as A and B. A 
particular case of a st ationary solutioii is the poloidal field 
configuration found bv lCumming et all (J2GG4), in which the 
electron velocity corresponds to a pure rigid-body rotation, 
so one has B — and an angular velocity = x^*-^ — 
constant everywhere. 

Aside from this very special purely poloidal ini- 

tial field {A 0,B = 0) will generate a toroidal compo- 
nent as well, because its electron velocity corresponds to 
a differential rotation that stretches the field lines in the 
azimuthal direction. However, the converse is not true: If 
the initial field is purely toroidal, i. e., initially A is ex- 
actly zero, it will remain t hat way, and only B will evolve 
(jUrpin fc Shalvbkovl . Il99lh . The evolution of such a field is 
much simpler than the general case, but still non-trivial, so 
we devote the next section to its study. 



4. Toroidal field 

4.1. Evolution: Burgers equation and current sheets 

The evolution of a purely toroidal magnetic field, now again 
including a resistive term with axisymmetric diffusivity, 
7y(i?, z), is given by 



dB 
'dt 



■ VB = R^V ■ (;^VB 



(10) 



wittQw = R^BVx'x'^4'- Thus, in the absence of dissipation 
(77 = 0), the quantity B can be viewed as being advected by 
the "velocity field" w like a scalar conserved quantity by 
a hydrodynamic fiow. This advection is very different from 
that of the magnetic flux B by the electron velocity field 
V — R^x^'t' X ^B, which can change the magnitude of the 
field by compressing or diluting the flux. The velocity held 
w is clearly perpendicular to B (i. e., poloidal and tangent 
to surfaces of x(-R, z) = constant). 

Figure 1 shows the surfaces of constant x for a particular 
neutron star model. The singular surface x = +00 is the 
union of the star's symmetry axis (on which i? = 0) and its 
surface (where n = 0). For all x larger than its minimum 
value Xmin , the "x-surfaces" are nested surfaces of toroidal 
topology, while Xmm defines an equatorial circle. 

In the absence of dissipation, the evolution of the mag- 
netic field at any given point depends only on its value at 
other points on the same X"Surface, not on the values on 
adjacent surfaces. We define a new coordinate s on each 
X-surface by the condition d/ds = R^Vx x V(/) ■ V, so a 
variation ds at constant x and (j) corresponds to a physical 
displacement \dr\ = R\Vx\ds. This allows the field evolu- 
tion to be written as the dissipationless Burgers equation 
(lBurgerslfl940LfT948h . 



Bt + BBs = 0, 



(11) 



^ There is an ambiguity in the choice of the vector field w, in 
the sense that adding to it another vector field that is everywhere 
perpendicular to VS does not change the evolution of 3. Our 
choice is convenient in the sense that the flow lines of w are 
independent of B and therefore do not change with time. 
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Fig. 1. A meridional cut of a spherical star with an assumed 
electron density profile n(r) = no[f — (r/ro)^], where r is 
the (spherical) radial coordinate, tq is the stellar radius, 
and no is the electron density at r = 0. Shown are "x- 
surfaces" corresponding to x/xo = 5, 20, and fO"*, where 
Xo = c/(47renorg). 

with the subscripts denoting partial derivatives. The do- 
main of this equation is a closed loop of constant x ^-nd (f>, 
therefore its bounda ry conditions must b e periodic. 

It is well known (|Burgersl Il940l |1948[ ) that, for a given 
initial condition S(s,0) = /(s), where / is an arbitrary 
function. Burgers' equation has an implicit, analytic solu- 
tion, 



(12) 



so the value B of the function at any point s is carried 
along the surface with "velocity" ds/dt = B. The larger 
values travel faster, so discontinuities are formed when they 
are about to "overtake" the more slowly moving, smaller 
values. 

The model considered here encompasses two previously 
studied models as limiting cases. One is a plane-parallel slab 
with a horizontal magneti c field and vertically de creasing 
electron density, in which IVainshtein et al.l ()200Cl[ ) showed 
that Burgers' equation governs the field evolution on hor- 
izontal planes. The other is a sphere or other axisymmet- 
ric body with a toroidal field and uniform e l ectron den- 
sity, in which simulations (jUrpin &: Shalvbkovl . Il991f ) yield 
that the magnetic field (initially taken to be symmetric 
with respect to the equatorial plane) progressively concen- 
trates in one hemisphere, where it develops a strong gradi- 
ent and eventually gets dissipated. This is easy to under- 
stand in terms of our results, as in this case the ^-surfaces 
are coaxial, cylindrical surfaces, along which the field drifts 
towards the surface, devel oping a discontinuity character- 
istic of Burgers' equation (|Araval l2002t iPons fc Geppertj . 
[200l . 

One feature of our, more general solution is not present 
in either of these special cases. Since the quantity B is con- 
served as it is carried along the x-surfaces, the magnetic 
field strength |B| = |SV0| ~ \B\/R changes in inverse pro- 
portion to the distance to the axis, R. In particular, when 
the flux near the surface of the star is carried from low lati- 
tudes to the vicinity o f a pole, the field strength increases, as 
found in simulations (jHollerbach fc Riidiged |2002| . |2004| ) . 



In both limiting cases mentioned above, one is forced 
to adopt somewhat arbitrary boundary conditions, in the 
first case because the slab is infinite, in the other because 
the electron density drops abruptly to zero at the surface 
of the star, where the x-surfaces end. This is not a problem 
in more realistic cases, in which the electron density decays 
continuously to zero at the stellar surface, and x-surfaces 
are closed and therefore have no boundaries. Of course, a 
real neutron star is not fully solid, but has a fluid core, 
which will be threaded by the x-surfaces. Thus, magnetic 
flux will generally be transported back and forth between 
the fluid core and the solid crust, an effect that will re- 
quire underst anding the Hall drift in flu i d mat ter, which is 
addressed bv iReisenegger fc Thompson! ()2007l ). 

In order to find the time of occurrence of the disconti- 
nuities in the magnetic field, one may consider the partial 
derivative of eq. p2)) with respect to s, Bg = {1 — Bst)f'{s~ 
Bt), from which one obtains Bg — /'/(I + tf). A disconti- 
nuity {Bs ±oo) is formed at time tdisc = l/rnax(— /'), 
at the comoving point where the initial spatial derivative 
had the largest, negative slope. (Note that the Hall effect 
is not invariant under spatial reflection.) 

These discontinuities in the magnetic field strength cor- 
respond to sufaces of infinite current density. In practice, 
before a "current sheet" becomes singular, resistive dissipa- 
tion must occur, even for arbitrarily small rj. This dissipa- 
tion occurs in a thin layer, whose width is ~ r]/{R\\/x\^B), 
where AB is the differen ce in the values o f B across the 
(near) discontinuity (Vain shtein et al.ll2000f ). However, the 
rate of dissipation depends exclusively on the rate at which 
the Hall effect transports the scalar variable B (related to 
the magnetic fiux) to the current sheet. 

4.2. Conserved quantities 

From the induction equation for a toroidal field (eq. [TU]). 
and using the identity 



V-(;^)=V.(VxxV0) = O, 



(13) 



we can show that any scalar function F{B) defines a "den- 
sity". 



_ F' 

p = w 

and a "flux" , 

that satisfy a continuity equation 



dp 
dt 



V- J 



i?2 



F"'{VBf 



(14) 



(15) 



(16) 



where primes (') denote derivatives of F with respect to its 
argument, B. Thus, the physical quantity whose density is 
p is transported by the Hall drift in the direction of w, i. e., 
along x-surfaces, and by Ohmic diffusion in the direction of 
decreasing B, while it is being destroyed by the right-hand- 
side term, which is also related to Ohmic diffusion. (For 
definiteness, in this discussion we assumed that B, F' , F" , 
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and F'" are all > 0.) Thus, in the absence of Ohmic diffu- 
sion {rj = 0), the volume integral over the toroid bounded 
by a x-surface. 



p{B[v,t]) dV, 



(17) 



x(r)<Xo 



is constant in time. We consider two specific examples of 
such conservation laws, namely those of magnetic flux and 
of magnetic energy. 

In order to obtain magnetic flux conservation, we must 
set F = B'^/2. For this choice, the right-hand-side of equa- 
tion (jl6p vanishes, indicating that magnetic flux can be 
transported by both Hall drift and Ohmic diffusion, but nei- 
ther created nor destroyed, and the continuity equation can 
be rewritten in terms of the magnetic flux density B = B/ R 
as 



5. Perturbations and instability 

5.1. Poloidal perturbations of the toroidal equilibrium field 

We now study the stability of this stationary, toroidal mag- 
netic fleld to a small, poloidal perturbation, Bi. The asso- 
ciated velocity fleld vi is toroidal, therefore vi x Bo = 0, 
and the linearized evolution equation for the perturbation 
reduces to 



dBi 



V X (vo X Bi), 



(23) 



dB 
'dt 



_d_ 

OR 
d_ 

dz 



1 



which implies that the fleld lines of the perturbation are 
carried along by the background electron flow fleld, vo = 
— (c/47rne)V x Bq, and the perturbation fleld Bi remains 
poloidal. 

Thus, it can be written as Bi = VA{x, s,t) x V0, with 
V^l <C \B\. Replacing into the above equation, one flnds 
that the potential evolves according to 



dB 
dz 



(18) 
(19) 



dA „ , dB dA 

— = -Vo • X-] — 

dt dx ds 



(24) 



i. e., the scalar quantity A is also carried along x-surfaces 
by the unperturbed electron flow. Thus, it can generally be 
written in terms of its initial condition as 



which can be interpreted as a two-dimensional continuity 
equation for B on the Cartesian R — z plane. 

For magnetic energy conservation, we set F — A{x, s,t) — A{x, s + x[dB / dx\t,Q). 
e3/(247r), obtaining 



(25) 



d_ 

dt 



' 1 fB\^' 






+ V- 


Stt \RJ 





w 



_nB_ 

47ri?2 



47ri?2 



{VBf. 



(20) 
(21) 



The right-hand-side is negative-definite, so magnetic energy 
can be destroyed by Ohmic diffusion, but not created. We 
note that the magnetic energy flux vector obtained here is 
generally different (in its Hall-drift part) from the standard 
Poynting vector. 



47r 



-E X B 



An 



VBxWcj) 



_qB_ 

47ri?2 



(22) 



However, their divergence, which determines their physical 
effects, is the same. 

In the scenario discussed in the previous section, in 
which Hall drift dominates and Ohmic dissipation is only 
signiflcant in thin current sheets, we expect no large-scale 
transport of magnetic flux or energy across x-surfaces, but 
only a redistribution on each surface by Hall drift and lo- 
cal Ohmic dissipation in the current sheets, which causes 
the magnetic energy to decrease while keeping the magnetic 
flux constant. Thus, the asymptotically resulting fleld con- 
figurations Bo = ;B(x)V(/) are minima of the magnetic en- 
ergy subject to conservation of flux in the region within each 
X-surface, and are the only toroidal equilibrium states as far 
as the Hall effect is concerned . We note that, like th e sta- 
tionary, poloidal fleld found bv lCumming et all ()2004l ). cor- 
responding to a rigidly rotating electron fluid, these conflg- 
urations are also not force-free, in the sense that j x B ^ 0. 
In the long term, of course, diffusion will allow magnetic 
flux transport across x-surfaces, eventually leading to mag- 
netic field decay. 



Initially, the smallness of Bi forces ^ to be a smooth 
function of position, taking similar values on adjacent 
points of different x-surfaces. However, when the circula- 
tion periods of the electrons around different x-surfaces are 
different, each point comoving with them will come close 
to others that were initially far away, so the values of A 
within a vicinity of fixed size around a given point will be- 
come progressively different, leading to a linearly increasing 
perturbation field. Another way of viewing this process is 
by realizing that the magnetic field lines of the perturba- 
tion are stretched as different parts are carried by electron 
currents circulating with different periods on different x~ 
surfaces. A simulation of this evolution is shown in Fig. 2. 

As the function A{Xj ^) is periodic, the perturbed 
field oscillates in time, periodically for the component per- 
pendicular to the x-surfaces, and with increasing ampli- 
tude for the tangential component. A similar, oscillat- 
ing behavior has been observ e d in n umerical simulations 
(Ho llerbach fc Riidieeil . l2002l 12004 IShalvbkov fc Urpinl . 
1996ir However, given their more complex geometries, a 
direct connection is difficult to establish. 



5.2. Plane-parallel analog 

In order to understand the relation of the instability 
found analytically by us to those found numerically by 
iRheinhardt fc GeppertI (|2002|) in a plane -parallel slab and 
later discussed by Gumming et al.l (j2004l ). we consider the 
plane-parallel geometry studied by these authors, in which 
a horizontal field varying nonlinearly with depth generates 
different velocities at different depths, stretching a small, 
vertical field component. 

We take the vertical coordinate to be z, and consider 
a background magnetic field Bo — f{z)x, which gener- 
ates an electron velocity Vq = —cf'{z)/[Aim{z)e]y, where 
primes denote derivatives with respect to z. If an additional, 
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Fig. 2. The evolution of a weak poloidal magnetic field 
component affected by a much stronger toroidal back- 
ground field. Shown are the poloidal field vectors on a cut 
through a star with electron density profile n(r) = rio[l — 
(r/rg)^], with central density rip = lO'^^cm"'^, and stellar 
radius tq ~ 10km. The background toroidal field has the 
form Bq = ;B(x)V(/), where we chose B{x) = boro{xo/x)^ ^ 
with bg = lO^^G, X = c/(47rnei?^), and xo its minimum 
value. The initial poloidal field was chosen uniform and 
pointing along the symmetry axis. The evolutionary times 
scale oc norQ/Bo- 



two sinusoidally varying potentials, as done in their equa- 
tions (4) and (5). For that representation, a magnetic field 
perturbation with vanishing horizontal wave vector would 
be identically zero. Nevertheless, in our solution as in that 
work, if the electron density is taken to be constant, the 
growth of the perturbation depends directly on the sec- 
ond spatial derivative of the background field strength, 
/"(z), suggesting that we are in fact seeing the same kind 
of instability. In our m odel, we confirm the conclusion of 
ICumming et al.l ((200l) that this finite second derivative is 
caused by a sheared electron flow velocity, which bends the 
vertical field component, creating an additional horizontal 
field. 

Curiously, in this simple geometry, we can actually cal- 
culate the nonlinear development of the instability. We con- 
sider a field with three components of arbitrary strength: 
Bx{z,t), which is initially the "background" field, but 
is now allowed to evolve arbitrarily large perturbations; 
the initially small (or vanishing), but growing component 
By{z,t)] and the vertical component Bz, which represents 
the initial perturbation. For this choice, the components of 
the induction equation become 



dB, 
dt 

dBy 
dt 

dBz 
dt 



d_ / cBz dBy 
dz \47rn(z)e dz 
d ( cBz dB, 



dz \47rn(z)e dz 



0. 



(26) 
(27) 
(28) 



These are the equations for "helicons" or "whistler waves" 
of arbitrary amplitude propagating along the z axis. For 
n(z) = constant, they have a dispersion relation oj = 
[cB z / {4,Trne)]k'^ , where to is the frequency and k is the mag- 
nitude of the wave vector (which points in the ±z direction). 
One particular solution is 



Bx{z,t) = BjnCOs{kz) cos{ujt), 
By{z,t) = B,nCOs{kz)sin{ujt), 
Bz ~ constant, 



(29) 
(30) 
(31) 



where the horizontal component initially points along the 
X axis, then turns to the y axis at a rate determined by Bz, 
and eventually turns around completely. Thus, the evolu- 
tion of the field in this case is strictly periodic (with the pe- 
riod equal to the Hall time scale of the "small" component 
Bz) and does not lead to dissipation on a Hall time scale. 
If Ohmic diffusion is introduced, the amplitude of the hori- 
zontal component will decay exponentially, on the resistive 
time scale l/{rjk'^). However, we suspect that this "clean" 
behavior is a peculiarity of this very symmetric configura- 
tion, and will not hold in more realisti c cases , including the 
simulations of iRheinhardt fc GeppertI (|2002f) . 



uniform, vertical field Biz is introduced, it causes the lin- 
ear growth of a component in the y direction, dBiy/dt = 
— \cBizf'{z)/{A'Kn{z)e)]'. In some sense, this is a partic- 
ular version or limiting case of the instabilities found by 
IRheinhardt fc GeppertI (|2002f ). with vanishing horizontal 
wave vector. However, it does not appear in their calcula- 
tion, as this field perturbation cannot be written in terms of 



5.3. Energy conservation in Hall instabilities 

In this section, we would like to clarify an issue that 
was rai sed, but in our view not full y clarified, by bot h 
Rhcinh ardt fc GeppertI (|2002f ) and Gu mming etall (|2004D . 
namely the conservation of energy in the growth of Hall 
instabilities. 
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In any slow evolution of a magnetic field, one can use 
the induction equation and Ampere's law to show that the 
magnetic energy density u = B"^ /{8tt) satisfies 



du 
'dt 



V-S 



(32) 



with the Poynting flux vector S — (c/47r)E x B. For pure 
Hall drift, the electric field, E = j x B/(ne), is perpendic- 
ular to the current, so the right-hand side vanishes and 
the magnetic energy satisfies a n exact conservation law 
(jGoldreich fc Reisenegg"eilll992f) . 

If we decompose the magnetic field, B = Bg-I-Bi, where 
Bo is the constant or slowly varying "background field" , 
whereas Bi is a small, more rapidly varying "perturbation" , 
each can be taken to satisfy an (exact) induction equation, 
d'Ba/dt — — cV X Eq, where 



Eo — jo X Bo, 

nec 

El = (jo X Bi + ji X Bo + ji X Bi) , 

nec 



(33) 
(34) 



and ja = (c/47r)V x Bq,, for a = 0, 1. The magnetic en- 
ergy density can be decomposed into three parts, u = 
uq + ui + U2, where ug = Bq/^Stt) and U2 = B\/{%'k) are 
positive-definite, while ui = Bq ■ Bi/(47r) is not. For uo, 
the derivation above can be retraced to show that it sat- 
isfies an exact conservation law. However, the evolution of 
ui and U2 yields 



dui 
'dt 
du2 
'dt 



V-Si 

V-S2 



I 



nec 



-jo • ji X Bi, 



-jo • ji X Bi, 



(35) 
(36) 



where Si = (c/47r)(Ei xBo+Eq xBi) and S2 = (c/47r)EiX 
Bi. Thus, the volume integrals of the two pieces, J ui 
and J U2 , are not individually conserved, but their sum, 
/ {ui 4-M2), is conserved. To linear order in Bi, only /ui is 
non-zero, and it is conserved. To quadratic order, J U2 can 
increase (in the case of an instability), being balanced by 
a corresponding decrease in /wi, which becomes progres- 
sively more negative. Thus, instabilities can only occur if, 
for some choices of the perturbation field Bi, the spatial 
integral of the right-hand side of equation ((36)) is positive. 
Clearly, this will not happen for current-free field configura- 
tions (jo = 0), for which this right-hand side vanishes iden- 
tically and no gr owth will occur. Somewhat more generally, 
ICumming et al.l |2004) have shown that for their stationary 
configurations, in which the current density is such that all 
the electrons rotate as a rigid body, the volume-integral of 
the right-hand side is zero for all perturbations that vanish 
on the boundaries, therefore also proving their stability. 

5.4. Current-free, poloidal background field 



An application of the result of § 15.31 is given by the follow- 
ing, plausible scenario for the evolution of the magnetic field 
in a neutron star. The field initially has an approximately 
axisymmetric configuration with both poloidal and toroidal 
components that are supported by currents that flow partly 
in the solid crust and partly in the fluid core of the star. In 
the core, the conductivity is very high, th e Hall drift is in- 
effective (jReisenegger fc ThompsonI 120071 ). and ambipolar 



diffusion not fast enough, so the fie ld is essentially frozen 
in some MHD-stable configuration (jBraithwaite fc Spruitl . 
2006). In the solid crust, Ohmic dissipation is more effec- 
tive, particularly if aided by Hall drift, and one might ex- 
pect the currents to decay, so the field might end up in a 
current-free configuration, which, by the arguments of the 
previous section, is stable under the Hall drift. 

In order to examine the evolution of the small pertur- 
bations, we use our earlier notation, with the background 
field written as Bo = V^o x Vi^, with Bq = (no toroidal 
component) and A*^o = (no current). It clearly does 
not evolve through Hall drift, since dAo/dt = 0. Adding a 
small perturbation, the latter will evolve as 



dBi 



1 d 



= i?^Bo • V(xAMi) = - — (xAMi), 
at X dQ 



(37) 
(38) 



where we have defined the coordinate C along field lines of 
Bo by the condition d/dC = xR^'Bq ■ V = (cBo/47me) • V, 
analogous to the definition of s on x-surfaces in § |4l 

One can combine the two equations, eliminating Bi and 
obtaining 



dt^ 



d_ 



ld_ 



(xAMi 



(39) 



This appears to be an hyperbolic differential equation with 
wavelike solutions corresponding to helicons or whistler 
waves travelling along the field lines of Bq. In the WKB ap- 
proximation (small- wavelength perturbations). Bo can be 
considered as uniform. A* reduces to the Laplacian, and 
this equation is exactly the whistler wave equation in a 
uniform background field. The behavior of long- wavelength 
perturbations may be more complex, but the arguments 
given above ensure that they will also be stable. 

The more gene r al, "u niformly rotating" configuration 
of ICumming et aTl (j2004l ) has Vq = n x r, and its field is 
therefore still poloidal, but not current- free. Since xA*^o — 
f2 = constant (although non-zero), its perturbations satisfy 
exactly the same equations ([571 through This is con- 
sistent with the energetic argument for it being stable as 
well. 

This means that the current-free poloidal field (and 
its generalization, the "uniformly rotating" field of 
ICumming et al.|[2004l ) is indeed stable under the Hall drift. 
The field in a neutron-star crust might be able to settle into 
such a state, which would remain in this form as long as 
the currents supporting it in the core do not change. 



6. Variational principle 

Since the Hall drift conserves magnetic energy, it is tempt- 
ing to search for a variational principle that would yield 
Hall equilibria, and in fact a fairly simple and natural one 
exists. Consider stationary points of the magnetic energy 
U — J B^/{8tt) subject to magnetic field perturbations 
i5B = V X (^ X B) that are due to an infinitesimal displace- 
ment field ^ that does not change the electron density, i. e., 
V • (n^) — 0. The latter condition implies that the displace- 
ment field can be written as ^ = (l/n)V x a, where a is 
an arbitrary vector field. This allows a perturbation of the 
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magnetic energy density to be written as 
B SB 



6u 



An 



-^-(1 X B) X B + — (j X B) X a 

An cn 



-a • V X 



j xB 



(40) 
(41) 

(42) 



Integrating over the volume, the divergence term becomes a 
surface integral that can be made vanish by requiring that 
the normal components of j, B and ^ all vanish on the sur- 
face. (The former two are required in order to ensure that 
the energy in the volume is conserved under Hall drift, the 
latter for consistency with the condition of not changing the 
electron density.) With these conditions, in order to have 
SU — for an otherwise arbitrary vector field a, we must 
have V X (j X B/n) =0, which is also the condition to have 
no Hall drift of the magnetic field. (It can be shown that, for 
these perturbations and the adopted boundary conditions, 
the magnetic helicity J A • B is automatically conserved, 
i. e., our constraint is stronger than that of requiring helic- 
ity conservation.) 

Put in a different way, our result states that, if the mag- 
netic field is such that it is not changed by the Hall drift, 
^ oc (l/n)V X B, then its energy will not change under the 
more general class of displacements £ — (l/n)V x a, with 
arbitrary a. 

Presumably, a subset of the configurations at which the 
magnetic energy is stationary, subject to the constraints 
discussed, will be maxima or minima of the magnetic en- 
ergy, i. e., to order or higher, all perturbations will cause 
magnetic energy changes of the same sign. For such config- 
urations, there will be only a small set of nearby configura- 
tions with a similar energy, and thus they will be stable un- 
der Hall drift. The other configurations at which the mag- 
netic energy is stationary will correspond to saddle points, 
inflection points, etc., whose energy is the same as that of 
a large set of other configurations that can be obtained by 
continuously deforming the magnetic field structure. We ex- 
pect that these will generally be Hall-unstable, in the sense 
that a small, initial perturbation can eventually make them 
evolve into a very different configuration. 



7. Conclusions 

We have found some interesting results regarding the Hall- 
drift evolution of an axially symmetric magnetic field in 
solid matter, in which electrons are the only mobile parti- 
cles: 

1) We have characterized "Hall equilibrium" configura- 
tions that do not evolve under the Hall drift (§[3]). 

2) We have found that an exactly toroidal field evolves 
discontinuities described by the Burgers equation, which 
dissipate on the Hall drift timcscalc fij 14.11). as i n the plane- 
parallel case studied earlier by Vainshte in et al.l (2000) . The 
field evolves into a stationary state that minimizes magnetic 
energy subject to flux conservation within a set of nested, 
toroidal "x-surfaces" \A.'2\i . 

3) However, we found these stationary, toroidal fields to 
be unstable to poloidal perturbations 15. ip . making them 
an unrealistic model for the fields of neutron stars. The 
instability involved is due to shearing of the perturbation 



neld by tne background ele ctron velocity, and tnus closely 
related to those studied by iRheinhardt fc GeppertI (|2002f ) 
and Gu mming et all (|2004l ). 

4) We have discussed the issue of energy conservation 
in Hall instabilities, giving a criterion for their occurrence 
and showing that current-free fields are stable 15. 3[) . We 
applied it to argue that a current-free, poloidal field in the 
crust, supported by axially symmetric currents flowing only 
in the core of the star, might represent a long-lived magnetic 
configuration of a neutron star 15. 4p . 

5) Finally, we showed that a Hall equilibrium is a sta- 
tionary point for the magnetic energy, subject to displace- 
ments that do not alter the electron density, and used this 
to argue that maxima or minima of the magnetic energy 
with respect to such perturbations will be stable under the 
Hall drift (§0. 

Acknowledgements. This work was supported by FONDECYT 
Regular Grants 1020840 (A. R. and J. P. P.), 1020844 (R. B. and 
P. A. A.), 1060644 (A. R.), and 1060651 (R. B.). J. P. P. is also sup- 
ported by a doctoral fellowship at Pontificia Universidad Catolica 
de Chile (PUC). A visit by D. L. to PUC was made possible by 
FONDECYT International Cooperation Grant 7020840. The authors 
thank R. Fernandez, U. Geppert, M. Lyutikov, H. Spruit, and C. 
Thompson for useful discussions. 



References 

Araya, P. A. 2002, Informe de Practica de Licenciatura (undergradu- 
ate thesis), Pontificia Universidad Catolica de Chile 
Biskamp, D., Schwarz, E., Zeiler, A., Celani, A., & Drake, J. F. 1999, 

Phys. Plasmas, 6, 751 
Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819 
Braithwaite, J., & Nordlund, A. 2006, A&A, 450, 1077 
Braithwaite, J., & Spruit, H. 2006, A&A, 450, 1097 
Burgers, J. M. 1940, Proc. Roy. Neth. Acad. Sci. (Amsterdam), 43, 2 
Burgers, J. M. 1948, Adv. Appl. Mech., 1, 171 
Gumming, A., Arras, P., & Zweibel, E. 2004, ApJ, 609, 999 
Geppert, U., & Rheinhardt, M. 2002, A&A, 392, 1015 
Geppert, U., Rheinhardt, M., & Gil, J. 2003, A&A, 412, L33 
Goldreich, P., & Reisenegger, A. 1992, ApJ, 395, 250 
HoUerbach, R., & Riidiger, G. 2002, MNRAS, 337, 216 
HoUerbach, R., & Riidiger, G. 2004, MNRAS, 347, 1273 
Jones, P. B. 1988, MNRAS, 233, 875 

Kulsrud, R. 2005, Plasma Physics for Astrophysics, Princeton: 

Princeton University Press 

Naito, T., & Kojima, Y. 1994, MNRAS, 2 66, 597 

Pons, J., & Geppert, U. 2007, preprint ( astro-ph/0703267 1 
Reisenegger, A., Prieto, J. P., Benguria, R., Lai, D., & Araya, P. A. 

2005, in Magnetic Fields in the Universe: From Laboratory and 

Stars to Primordial Structures, AIP Conference Proceedings, vol. 

784, eds. E. M. de Gouveia dal Pino, G. Lugones, & A. Lazarian, 

p. 263 

Reisenegger, A., & Thompson, C. 2007, in preparation 
Rheinhardt, M., & Geppert, U. 2002, Phys. Rev. Lett., 88, 101103 
Rheinhardt, M., Konenkov, D., & Geppert, U. 2004, A&A, 420, L33 
Shalybkov, D. A., & Urpin, V. A. 1996, A&A, 321, 685 
Urpin, V. A., & Shalybkov, D. A. 1991, Sov. Phys. JETP, 73, 703 
Vainshtein, S. L, Chitre, S. M., & Olinto, A. V. 2000, Phys. Rev. E, 
61, 4422 



